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RENORMALIZED SELF-INTERSECTION LOCAL TIME FOR 
FRACTIONAL BROWNIAN MOTION 

By Yaozhong Hu 1 and David Nualart 2 

University of Kansas and Universitat de Barcelona 

Let Bf be a d- dimensional fractional Brownian motion with 
Hurst parameter H £ (0, 1). Assume d > 2. We prove that the renor- 
malized self-intersection local time 

1=1 ( 5{B" -B")dsdt-m( f ( S(B" - B")dsdt\ 
Jo Jo vo Jo ' 

exists in L 2 if and only if H < 3/ (2d), which generalizes the Varadhan 
renormalization theorem to any dimension and with any Hurst pa- 
rameter. Motivated by a result of Yor, we show that in the case 3/4 > 
H > , r(e)£ s converges in distribution to a normal law N(0, Ta 2 ), as 
e tends to zero, where £ e is an approximation of £, defined through (2), 
and r(s) = | logep 1 if H = 3/ (2d), and r(e) = e d ~ 3/(2H) if 3/(2d) < 
H. 



1. Introduction. The fractional Brownian motion on R rf with Hurst pa- 
rameter H G (0,1) is a d-dimensional Gaussian process B H = {B^,t > 0} 
with zero mean and covariance function given by 

E(B^B^) = 5 -f{t 2H + s 2H -\t- s\ 2H ), 

where i,j = l,...,d, and s,t > 0. We will assume that d>2. The self- 
intersection local time of B H is formally defined as 

(1) /= [ T C S (B t H - B^)dsdt, 

Jo Jo 
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where §o(x) is the Dirac delta function. It measures the amount of time 
that the process spends intersecting itself on the time interval [0,T] and 
has been an important topic of the theory of stochastic process. A rigorous 
definition of this random variable may be obtained by approximating the 
Dirac function by the heat kernel 

p £ (x) = (27rey d / 2 e W {-\x\ 2 /(2e)}, 

as e > tends to zero. We denote the approximated self-intersection local 
time by 

(2) I £ = f T fp £ {B?-B^)d S dt, 

Jo Jo 

and a natural question is to study the behavior of I £ as e tends to zero. 

For H = 1/2, the process B H is a classical Brownian motion. The self- 
intersection local time of the Brownian motion has been studied by many 
authors [see Albeverio, Hu and Zhou (1997), Calais and Yor (1987), He, 
Yang, Yao and Wang (1995), Hu (1996), Imkeller, Perez-Abreu and Vives 
(1995), Varadhan (1969), Yor (1985) and the references therein]. In the case 
of the planar Brownian motion, Varadhan (1969) has proved that I £ does not 
converge in I? but it can be renormalized so that I £ — log(l/e) converges 
in I? as e tends to zero. The limit is called the renormalized self-intersection 
local time of the planar Brownian motion. This result has been extended by 
Rosen (1987) to the (planar) fractional Brownian motion, where it is proved 
that for 1/2 < H < 3/4, I £ - C H Te^ l+l / {2H "i converges in L 2 as e tends to 
zero, where Ch is a constant depending only on H [see also (4): Ch = Ch,2\- 

On the other hand, for classical Brownian motion and when the dimension 
d > 3, Yor (1985) and Calais and Yor (1987) have proved that the random 
variables 

r(log(l/e))- 1 /2(j e _E(J £ )), if d = 3, 

\ £ <V2-3/2 (/£ _ E(4))) if d>3j 

converge in law to normal distribution (with a different approximation). 

The first aim of this paper is to extend Rosen's and Varadhan's results 
to arbitrary dimensions and with Hurst parameter H < 3/4. Along this line, 
we have established the following result. 

Theorem 1. Let I £ be the random variable defined in (2). We have: 

(i) If H <l/d, then I £ converges in L 2 as e tends to zero. 

(ii) Ifl/d<H<3/(2d), then 

(3) I £ -TC H4 E- d ' 2+l IW 
converges in L 2 as e tends to zero, where 

(4) C H4 = (2n)- d / 2 (z 2H + l)- d / 2 dz. 

Jo 
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(iii) Ifl/d = H<3/(2d), then 

(5) h T w/o log(l/e) 

V ; 2H(2^) d / 2 6V 1 ' 

converges in L 2 as e tends to zero. 

(iv) If H > 3/ (2d), then the difference I £ — E(/ e ) does not converge in I? . 

That means, if H < 3/ (2d), the difference I £ — E(J e ) converges in L 2 as e 
tends to zero to the renormalized self-intersection local time. 

In the case H > 3/ (2d) we have established the following version of the 
central limit theorem. 



Theorem 2. Suppose 3/ (2d) < H < 3/4. Then the random variables 

( , f (logCl/e))" 1 /^ - E(J £ )), i/ff = 3/(2d), 

1 j U d / 2 - 3 /(^)(/ £ -E(/ £ )), ifH>3/(2d), 

converge as e tends to zero in distribution to a normal law N(0,To~ 2 ), where 
a 2 is a constant depending on d and H . 

We conjecture that, as functions of the terminal time T, these processes 
converge in law to {o~Bt,T > 0}, where Bt is a Brownian motion. We plan 
to discuss this problem in a forthcoming paper. 

For H > 3/4 our method of proof of Theorem 2 does not work and the 
study of the fluctuations of I £ — E(/ e ) in this case is an open problem. 

The proof of Theorem 1 is based on some estimates deduced by Hu in 
(2001) from the local nondeterminism property of the fractional Brownian 
motion. The proof of Theorem 2 is more involved, and the main ingredient 
is to show that the projections on each Wiener chaos converge in law to 
independent Gaussian random variables. The proof of these convergences 
is based on a recent general criterion for the convergence in distribution 
to a normal law of a sequence of multiple stochastic integrals established 
by Nualart and Peccati (2005) and by Peccati and Tudor (2003), in the 
multidimensional case. We have extended their results to a slightly more 
general setting applicable to the renormalization of self-intersection local 
time of fractional Brownian motion and to this end we have established a 
general central limit theorem which has its own interest. 

In the case of the classical Brownian motion, the convergence in law of 
the Wiener chaos projections of the random variables appearing in (6) has 
been proved in de Faria, Drumond and Streit (2000) using a martingale 
approach which cannot be extended to the case H ^ 1/2. The corresponding 
convergence in law to a family of independent Brownian motions has been 
established in Rezgui and Streit (2002). 
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The paper is organized as follows. In Section 2 we prove Theorem 1. Sec- 
tion 3 is devoted to a general central limit theorem using chaos expansions. 
Some examples are given and will be used to deal with the self-intersection 
local time of the fractional Brownian motion. Section 4 deduces the chaos 
expansion of the approximated self-intersection local time, and in Section 5 
we prove Theorem 2. The proofs of the main results are based on a sequence 
of technical lemmas, stated and proved in Section 6. 

2. Renormalized self-intersection local time, case dH < 3/2. Let B H = 

{Bt ,t > 0} be a ci-dimensional fractional Brownian motion of Hurst param- 
eter H £ (0,1). 

Consider the approximation I e of the self-intersection local time intro- 
duced in (2). From the equality 

p £ (x) = (27r)- d ^exp{i(e,x)}exp|-^|^ 
and the definition of I e , we obtain 

(7) I £ = (2vr)- d f T f ( exp{i(Z,B?-B?)}exJ- £ -^-XdZdsdt. 

Jo Jo JR d I 2 J 

Therefore, 

K(I 2 ) = (2vr)~ M / / E(exp{i£(Bf - Bf )} exp{«/(i# - B*)}) 

(8) JTk2d 

x exp{-(e|£| 2 + e\r]\ 2 )/2} d£ drj ds dt ds' dt', 

where 

(9) T = {{s, t, s', t') : < s < t < T, < s' < t' < T}. 

Throughout this paper we will make use of the following notation, for any 
T=(s,t,s',t>): 

(10) \( T ) = \t-s\ 2H , P (r) = \t' -s'\ 2H 
and 

(11) M (r) = \[\s- t'\ 2H + |a' - t\ 2H -\t- t'\ 2H -\s- s'\ 2H \. 

Notice that A is the variance of B^ 1,1 — B^' 1 , p is the variance of B^' 1 — B^' 1 

H 1 if 1 H 1 H 1 

and n is the covariance between B t ' — By' and B t , ' — B ,' , where 
B H > 1 denotes a one-dimensional fractional Brownian motion with Hurst pa- 
rameter H. 

With this notation, for any £,rj € M. d , we can write 

(12) E[«£,£f - Bf) + { V ,B$ - B%)) 2 ] = \\tf + p \rf + 2 p{t, v ). 
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As a consequence, from (8) and (12) we deduce for all e > 0, 

E(if) = (27r)- 2d [ / e -((^)l«l 2 + 2 ^)+(P+-)l^l 2 )/2 ded??dr 

(13) JtUm 

On the other hand, the expectation of the random variable I e is given by 

E(J e )= [ T f p £Mt _ s?H {Q)dsdt 
Jo Jo 

(14) =(2vr)- d / 2 [ T [\e + \t-s\ 2H )- d ' 2 dsdt 

Jo Jo 

= (2vr)- d / 2 [ T (T-s)(e + s 2H )- d / 2 ds. 
Jo 

PROOF of Theorem 1. Assertion (i) follows easily from (13) and (14). 
From (14), making the change of variables s = ze l ^ 2H \ we obtain, if 1/d < 
H<3/(2d), 

-l/(2H)-d/2 rTe 1 /^ 

EC) = "-^ym- 1 (t - „v<w )(j » + ir ,/ 2 dz 
=E 1 '«-« ! re w+ »(e). 

For H = 1/d we get 

rio g (i/g) , 

E(i e ) = ; 73777 + o(e). 

Hence, the convergence in L 2 of the random variables (3) and (5) is equiva- 
lent to the convergence of I e — E(J e ). 
From (13) and (14) we obtain 

E(I £ I V ) - E(J e )E(J„) 

= (27r)- d / [((A + e)(p + v)- fT d ' 2 -((\ + e){p + r?))"^ 2 ] dr. 



IT 

Therefore, a necessary and sufficient condition for the convergence in L 2 of 
I e -E(J e ) is that 

(15) ~ r =: jfj(Ap - /x 2 )"^ 2 - (Ap)^/ 2 ] dr < oo. 

Finally, Lemma 11 allows us to conclude the proof of the theorem. □ 
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3. Central limit theorem via chaos expansion. In this section we will es- 
tablish a general central limit theorem for nonlinear functionals of a Gaussian 
process, based on the Wiener chaos expansion. The proof of this theorem 
uses some recent results by Nualart and Peccati (2005) and by Peccati and 
Tudor (2003). 

Consider a separable Hilbert space H. Let {ek,k > 1} be a complete or- 
thonormal system in H. For every n > 1, we denote by H &n the nth sym- 
metric tensor product of H. For every p = 0, . . . ,n and for every / G H &n , we 
define the contraction of / of order p to be the element of ^® 2 ( n_ p) defined 
by 

oo 

f®pf= ^' e *i ® •••® e ; P )tf®p ® (/,e;i ® •■•® 

il,...,ip=l 

Suppose that A" = {X(/t), /i G if} is an isonormal Gaussian process on H. 
This means that X is a centered Gaussian family indexed by the elements 
of H, defined on some probability space (0,.F, P) and such that, for every 
h, ti £ H, 

E(X{h)X{ti)) = (h,ti) H . 

Assume that T is a c-algebra generated by X. 

For every n > 1 we will denote by I n the multiple Wiener integral [Hu 
and Kallianpur (1998)]. 

Theorem 3. Consider a sequence of square integrable and centered ran- 
dom variables {Fk,k > 1} with Wiener chaos expansions 

oo 

Fk = X! In(fn,k)- 
n=l 

Suppose that: 

(i) lunjv^oolimsup^^E^Jv+i^'ll/n^ll^n = 0; 

(ii) for every n>\, lim fc ^oo rc!||/n,fc||#®n =o- 2 n ; 

(iii) E^=i^ = ^ 2 <oo; 

(iv) for all n> 2, p= l,...,n- 1, lim^oo ||/ n;fe ® p f n ,k\\ 2 H mn-p) =°- 

Then converges in distribution to the normal law N(0,a 2 ) as k tends 
to infinity. 

Proof. By Theorem 1 in Nualart and Peccati (2005), conditions (ii) and (iv) 
imply that for each fixed n > 1 the sequence of multiple integrals I n (fnk) 
converges in distribution to the normal law iV(0,<7^) as k tends to infinity. 
Furthermore, by Theorem 1 in Peccati and Tudor (2003), for each n > 1, 
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the vector (h(fi : k), ■ ■ -,In{fn,k)) converges in law to the vector (£ l5 . . . ,£ n ), 
where {£, n , n > 1} are independent centered Gaussian random variables with 
variances {cr^,n > 1}. Now let / be a C 1 function such that |/| and \ f'\ are 
bounded by 1 . For every TV > 1 , set 

N 
n=l 

Set C (N) = En=i £», and £ = J2n=i £«■ Then 
|E(/(F fc ))-E(/(e))| 

<|E(/(F fc ))-E(/(Ff } ))| 

+ |E(/(Ff } )) - E(/tfW))| + |E(/(^)) - E(/(0)| 

/ oo \ 1/2 

< E ^ll/Mlllr^ +|E(/(Ff ) ))-E(/(^))| 

\n=7V+l / 

+ |E(/(^))-E(/(0)|- 

Then, taking first the limit as k tends to infinity, and then the limit as N 
tends to infinity, and applying conditions (i), (iii) and the convergence in 
law of the vector (^(/i^), . . . ,I n (fn,k)), we get the desired result. □ 

Remark 1 . A sufficient condition for (i) and (iii) is that for all n > 1 , 
and k > 1 , 

n !||/n,fc||/f®n < S n , 

where E^i^n < oo. 

Remark 2 . The assumptions of Theorem 3 imply the convergence in law 
of the projections on the Wiener chaos to a family of independent Gaussian 
random variables. 

The following propositions contain examples of applications of the above 
central limit theorem which will be useful in proving the renormalization 
result for the fractional Brownian motion when dH > 3/2. 

Let TC be the Hilbert space defined as the closure of set 8 of step functions 
on R + with respect to the scalar product 

{ 1 [0,t], 1 [0,s])H = 2 + s )■ 

Let B H = {B^ ,t > 0} be a (2m)-dimensional fractional Brownian motion 
with Hurst parameter H, with m > 1. Then, B H is an isonormal Gaussian 



8 



Y. HU AND D. NUALART 



process associated with the Hilbert space H = TL . We denote by /„ the 
multiple stochastic integral which provides an isometry between the sym- 
metric tensor product (7i 2m ) & equipped with the norm \/n7|| ■ ||#®n and the 
nth Wiener chaos of B H . 

We will make use of the following notation: 

Ki{x,y,z) = \z + y\ 2H + \z - x\ 2H - \z + y-x\ 2H - z 2H , 

K 2 (x,z) = Ki(x,x,z). 



Proposition 4. Let {Bi^ , t > 0} be a (2m) -dimensional fractional Brow- 
nian motion with Hurst parameter H G (0, 3/4). For any < s < t, we set 

2m 

(i6) z Sti Hi />•"•' />',"•'). 

i=l 

Define, for any fixed x > 0, 

i-T 



1 f T 

At{x) = —== Z s , s+X ds. 
VT Jo 



Then as T tends to infinity At(x) converges in distribution to a normal 
law iV(0, G(x, x)), where 



poo 

G(x,x)=2 1 - 2m K 2m (x,y)dy. 
Jo 



Proof. We have 

2l-2m pT pt 



ni — zm pi pt 

E[A 2 T (x)} = —^ / K 2m {x,t-s)dsdt 
J Jo Jo 

= 2 l-2m j T K 2m { x,y){l-y/T)dy 
Jo 



and, clearly, this converges to G(x,x) as T tends to infinity. We can write 

A T (x) =hm{f2m > T), 

where 

-y pX 2m 

f2m,T(hm,ri,...,r 2m ) = ——— 7 = / ~[l[ s ,s+x]( r i)ds, 

{2m)WT Jo i=1 

if all components of the multi-index i2 m = (ii, ■ ■ ■ , i2m) axe different and zero 
otherwise. Then, by Theorem 1 of Nualart and Peccati (2005), it suffices to 
show that for any p, 1 < p < 2m — 1, 

(17) T^oo ^ 2m ' T ® p /2m,T||^®(2m- P ) = 0. 
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The contraction of order p of the kernel f2m,T is given by 

(/2m,T ®p f2m,T){hm-p, H., . . . , r 2m _ p ,hm-p, U\ , . . . , U 2m - p ) 
, ji T 2m-p 

= (( 2m )!)2 r2 p X y n w^^w^)*^*- 

In order to show (17) it suffices to prove that 

lim - / K$(x,h - h)Kl(xM ~ h) 

T->CO 1 J0<tl<t2<t3<*4<T 

(18) 

x K 2m ~ p (x, t 3 - t 2 )K 2m - p (x, t 4 - h) dt = 0. 
By L'Hdpital's rule, (18) holds if 

lim i / iff (x, t 2 - h)K$(x, T - t 3 ) 

(19) x K% m ~ p (x,t 3 -t2)K% m ~ p (x,T -t 1 )dt 1 dt 2 dt 3 
= 0. 

Set t 2 — t\ = a, t 3 — t 2 = b, T — t 3 = c. In order to show (19) it suffices to 
show that 

lim ± / K p (x, a)K p (x, c)K 2m ~ p (x, b) 

(20) 2m 

x K 2 m p (x,a + b + c)dadbdc = 0. 

To show (20) we will make use of the following inequalities: 

(21) \K 2 (x,a)\ < k H a 2H ~ 2 x 2 , 
if x < § , where k H = 2 3 ~ 2H H(2H - 1), 

(22) a + b + c>3(abc) 1/3 
and 

(23) \K 2 (x,a)\ < 2{ax) H . 

We decompose the set G = {0 < a + b + c < T} as follows: 

8 

G=\jG h 

i=l 
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where 



Gi 


= { 


a 


b,c) 


eG 


a > 2x, b > 2x, 


c > 2x} 


G 2 


= { 


a 


b,c) 


eG 


a > 2x, b > 2x, 


c < 2x} 


G 3 


= { 


[a 


b,c) 


eG 


a > 2x, b < 2x, 


c > 2x} 


G 4 


= { 


[a 


b,c) 


eG 


a < 2x, b > 2x, 


c > 2x} 


G 5 


= { 


[a 


b,c) 


eG 


a > 2x, b < 2x, 


c < 2x} 


G 6 


= { 


a 


b,c) 


eG 


a < 2x, b > 2x, 


c < 2x} 


G 7 


= { 


a 


b,c) 


eG 


a < 2x, b < 2x, 


c > 2x} 


Gs 


= { 


a 


b, c) 


eG 


a < 2x, b < 2x, 


c < 2x} 



Set, for i = 1, . . . , 8, 

Air = ^ [ Kl{x,a)KUx,c)Kl m - p {x,b)Kl m - p {x,a + b + c)dadbdc. 

In what follows, we will denote by A; a generic constant that may depend on 
H and x. 

Case 1. Consider first the integral over the region Gi. Then, using (21) 
and (22) we obtain, for T large enough, 

K\<- [ a (^-2) Pb (2H-2) Pc (2H-2)(2m- P ) {a + & + c) (2ff-2)(2m-p) 

T Jgi 

< t. (£ a 4/3(2H-2) P d \ 2 a 4/3(2H-2)(2 m - P ) ^ 

which converges to zero as T tends to infinity, because 

j(2H - 2)(2m +p) + 2 < |(2i? - 2)3 + 2 = 8if - 6 < 0. 

Case 2. Using (21), (22) and (23) we obtain 

A 2 T<^ I a^ 2H - 2 H^ 2H - 2 ^c H ^ m - p \a + b + c)^ H - 2 ^ 2m - p Uadbdc 
T Jg 2 

< kT^ H ~ 4)i - m+p/2)+1 , 
which converges to zero as T tends to infinity, because 

{AH - 4)(m + p/2) + 1 < 3(2if - 2) + 1 = 6H - 5 < 0. 
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Case 3. Using (21), (22) and (23) we obtain 

A 3 T <t f a (2H-2) PbHPc (2 H -2)(2 m - P ), +b + c) (2H-2)(2rn- P ) dadbdc 
< ^^(4^-4)(2m-p/2)+l o 

because 2m — p/2 > 2m — (2m — l)/2 = m + 1/2 > 1, and, so 
{AH - 4)(2m - p/2) + 1 < 4H - 3 < 0. 

Case 4. The integral over the set G4 is handled as in Case 3. 

Case 5. Using (21), (22) and (23) we obtain 

A 5 <A / a ^-2)p b H Pc H(2m- P ) ia + b + c) (2H-2) { 2 m - P ) dadbdc 

T Jg 5 
= kT^ H -^ m 0. 

Finally, the integrals over the sets Gg, G7 and Gs are treated in a similar 
way. 

□ 

Proposition 5. Consider the processes A^(x) defined in Proposition 4. 
If3/(2d) < H < 3/4, i/ien t/ie random variables 

roc 

Y T := / (l + x 2H )- d / 2 - m A T (x)dx 
Jo 

converge in law to the normal distribution N(0,a 2 ), as T tends to infinity, 
where 

roc roo 

(24) a 2 := / / (1 + x 2H y d l 2 - m (\ + y 2H )- d / 2 - m G(x, y) dxdy 
Jo Jo 

and 

G(x,y):=2 1 - 2m K 2m (x,y, z) dz. 
Jo 

Proof. Notice that 

E[A T (x)A T (y)] = —^ / K 2m (x,y,t-s)dsdt 
1 Jo Jo 

= 2 l-2m [ T K f™ {x ,y jZ){ i- z/T)dz , 
Jo 

and this clearly converges as T tends to infinity to G(x,y). 
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As a consequence, Proposition 4 and Proposition 2 of Peccati and Tu- 
dor (2003) imply that the finite-dimensional distributions of the process 
{At(x),x > 0} converge to the finite-dimensional distributions of a Gaus- 
sian centered process with covariance function G(x,y). In order to complete 
the proof of the lemma we need to show the following facts: 

(i) 

poo 

(25) / (l + x 2H )- a! / 2 - m supE(|A T (x)|) <oo. 

JO T 

(ii) For all K > 0, 

(26) lim E[(A T (x) - A T (y)) 2 ] -> 

T — >oo 

as \x — y\ < 5 — > 0, and x, y < K. 
(hi) a 2 < oo. 

Proof of (i). Property (25) follows from 

r^ + x^y^-^Gix^dx 

Jo 

poo / p2x \ 1/2 

<J o {l+x 2H r d ' 2 ~ m ^J^ K 2m (x,y)dyj dx 

poo / poo \ 1/2 

(27) + j (l + x 2H )- d / 2 - m ^ K 2m (x,y)dyj dx 

< k H , m jH(l + x ^ rd ,2- m ^J 2 J {y AH m + x W m) ^ V2 ^ 

+ fe,m / o °°(l + x 2 ")^/ 2 "™ QT y( 2H ~V 2m x* m dy^j V2 

;>oo 

<k' Hm (1 + x 2H)-<i/2-m x 2mtf + l/2 dx < ^ 

' Jo 

because Hd > 3/2. □ 

Proof of (ii). We have 
supE[(A T (x)-A T (y)) 2 ] 

T 

poo 

<2 l ~ 2m / \K 2m (x,z) + K 2m (y,z) - 2K 2m (x,y,z)\dz, 
Jo 

and, by dominated convergence this tends to zero as \x — y\ tends to zero in 
the compact interval [0,if]. □ 



FRACTIONAL SELF-INTERSECTION LOCAL TIME 13 

Proof of (iii). This follows from (27) and the fact that, by Cauchy- 
Schwarz inequality, we have G(x,y) < y / G(x,x)G(y,y). □ 

Let us now show that (i)-(iii) imply the convergence in law of Yt to 
£, where £ is a random variable with the iV"(0,<7 2 ) distribution. Let / be a 
function such that |/| and \ f'\ are bounded by 1. We make the decomposition 

\E(f(Y T ))-E(f(0)\ 

< \E(f(Y T )) -E(f(Y± K) ))\ + |E(/(y^)) -E(/(^)))| 

+ | E (/(^)))-E(/(0)|, 

where 

4 K) = / (1 + x 2H y d / 2 ~ m A T (x) dx, 
Jo 

^ has the law N(0,a 2 K ), and 

a 2 K = 2 l-2m f K f K ^ + x 2H ) -d/2- m{l + y 2H ) -d/2- mG ^ y)dxdy _ 

Jo Jo 

Property (i) implies that lim^oo sup T \E(f(Y T )) - E(f(Y^ K) ))\ = 0. 

Property (iii) implies that lim^oo \E(f(^)) - E(/(f))| =0, and prop- 
erty (ii) implies that for any fixed K> 0, lim T ^oo \E(f(Y^ K) ))-E(f(^ K> >))\ = 
0. 

This completes the proof of the proposition. □ 

Proposition 6. Assume H = 3/ (2d) < 3/4. Consider the processes 

X T = -^J_ jT(l + x 2 ^)- d / 2 — (j^ Z SiS+x (fa) dx, 

where T > 0, and Z SjS + x /ias 6een defined in (16). T/ien converges in 
distribution, as T tends to infinity, to a normal law N(0,o~ 2 ), where 

poo poo 

a 2. = 2 l-2 m / ( xy )-Z/2-2rnH K 2, m ^ yA)dxdy 

Jo Jo 

Proof. As in the proof of Proposition 4 we can show that, for each 
fixed x > 0, the random variables , 1 iT Z* ,+ T ds converge to in 1? 

' yj T log T J0 s ' s+x 6 

as T tends to infinity. For this reason we cannot apply the same method of 
proof as in Proposition 4, and new ideas are required. The basic ingredient of 
the proof will be a scaling argument and Theorem 1 of Nualart and Peccati 
(2005). 
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Let us see first that the asymptotic behavior of Xt as T tends to infinity 
is the same as Yp, where 

yt ~ tt^t r x ~ 3 ' 2 ' 2mH {[ z >,°+* d, >) ^ 

In fact, 

n\x T -y t \]< jT°°(i +x 2H r 5 / 2 -™E(\ l £ z s , s+x d S ^ dx, 

and using the estimate (27) it is not difficult to see that this converges to 
zero as T tends to infinity. 

So, it suffices to show that Yt converges in distribution as T tends to 
infinity to the normal law iV(0, a 2 ). The proof will be done in several steps. 

Step 1. We claim that 
(28) lim E(y|) =a 2 < oo, 

T — >oo 

where a 2 is defined in (24). In fact, 

9 poo poo pT pt 

E(Y 2 ) = —— / (xyy 3 / 2 - 2mH n(s,s + x,t,t + y)d S dtdxdy 

TlogT J i J i Jo Jo 



2 



1— 2m poo poo pT / 7 \ 

J i J (xy)- 3/2 - 2mH J o K 2m (x,y,z)(l--)dzdxdy. 



logT 

Thus, by L'Hopital's rule 



poo poo 

lim E(Y$) = lim T2 1 ~ 2m / / (xy)~ 3 ^ 2 ~ 2mH Kf m (x, y, T) dx dy 

^oo T^oo Ji Ji 

poo poo 

lim 2 l ~ 2m / / (xy)- 3/2 - 2mH K 2m (x,y,l)dxdy 

"-^oo Jl/TJl/T 



T- 

a 2 . 



Step 2. Consider the decomposition 
(29) Y T = Y TA + Y Tt2 , 

where 

r T/2 
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From the scaling and the stationary increments properties of fractional 
Brownian motion, it follows that the random variables Yt \ and Yt,2 have 
the same distribution as Yrs, where 

Y T , 3 = —± = r x -^au( f T Z ds ) dx . 



v / 2TlogT J 2 

We are going to compute the second and fourth moments of Yt, using the 
decomposition (29). We have 

E(y£) = E(y| >1 ) + E(y| j2 ) + 2E(y Tjl y Ti2 ) 
= 2E(y£ 3 ) + 2E(y T)1 y Ti2 ). 

Clearly, limT^ 00 E(y^ 3 ) = a 2 /2, and this implies that 

(30) lim E(y Ti iy Ti2 ) = 0. 

T — >oc 

Consider now the fourth moment of Y?: 

(31) e(y£) = 2E(y£ 3 ) + 4E(y^ 1 y T , 2 ) + 4m.{y t ^y^ 2 ) + 6E(y| )1 y^ i2 ). 

Step 3. We claim that 

lim E(y T 2 1 y T 2 2 ) = a 4 /4, 

(32) 

lim E(y^ 1 y Ti2 ) = 0, lim E(y T> iy^ 2 ) = 0. 

Let us prove the first identity of (32). The random variables (ir,i^T,2) have 
the same joint distribution as (Zt,iZt,2), where 

z Til = ~=l= r x -w-*™h( f i2 z ds \ fa 

VlogTJi/T \Jo J 

z Ti2 = -±- r x -m-^n x z ds ) dx , 

VlogTJl/T \Jl/2 J 

As a consequence, we have 

(33) E(y T 2 ;i y T 2 2 ) = (logT)- 2 ^E^_Q A^ dfi, 

where Aj = - B tl , for i = 1, 2, 3,4,2?= [1/T, oo) 4 x [0, 1/2] 2 x [1/2, l] 2 

and d/j, = Ylt=i x i 3 ^ 2 2mH dxdt. Using the notation dij = E(AjAj) for i, 
j = 1, ... ,4, we can write 

/ 4 \ 



E II Al = dl2d 34 + di 3 d 2 4 + d u d 2 3 
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and 

(34) Wll^ll = (di2d 34 ) 2m + R- 



We have 



(logT)" 2 f (d 12 d 34 ) 2m dfi = E(Y 2 ^)E(Y 2 ^, 
Jv 



and this converges to <r 4 /4 as T tends to infinity. On the other hand, the 
residual term R in (34) does not contribute to the limit (33). In fact, by 
Holder's inequality we obtain, for any indexes (i,j,k) such that i + j + k = 
2m, and j > 1 or k > 1, 



(logT) 



(C?12C?34) 1 (^13^24^ (^14^23) 4* 

<(logT)- 2 ^(d 12 d 3 4) 2m V 

/• „ V'/(2m)/ , xfc/(2m) 

J v (d 13 d 24 ) 2m d t x\ (j^d u d 23 ) 2m dv) 



Then, this converges to zero as T tends to infinity because of (30), j + k > 1 
and 

(logT)- 2 / (d 14 d 23 ) 2m d^ = (logT)" 2 f (d 13 d 24 ) 2m d^ = (E(Y Tjl Y T . 2 )) 2 . 
Jv Jv 

This completes the proof of the first identity of (32). The proof of the other 
two identities of (32) follows the same argument. 

Step 4. By the same arguments as above and Minkowski inequality it 
follows that 

E(Y4) = (logT)- 2 J^(d 12 d u + d l3 d 2A + d 14 d 23 ) 2m dfi 

<3 2m (logTr 2 J £ (d 12 d 34 ) 2m dv 
= 3 2m (E(Y^)) 2 , 

where £ = [1/T, oo) 4 x [0, l] 4 . So, the sequences Yp and Yt,3 are bounded 
in L 4 , and 

(35) E((Y T - Y Tj3 ) 4 ) < 3 2m (E(Y T - Y T , 3 ) 2 ) 2 , 

which clearly converges to zero as T tends to infinity. Then, from (35), (31) 
and (32) we obtain lim^oo E(Y^) = 3cr 4 . Finally, the result follows from 
Theorem 1 of Nualart and Peccati (2005). □ 
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4. Wiener chaos expansion of the self-intersection local time. In this sec- 
tion we will first compute the Wiener chaos expansion of the approximation 
of the self- intersection local time I e defined by (2). 

Given a multi-index i n = (ii, . . . ,i n ), 1 < ij < d, we set 

a(i n )=E[X H ---X in ], 

where the Xi are independent iV(0, 1) random variables. Notice that 

^_ (2mi)!---(2m d )! 
aM_ N!-N!2-' 

if n = 2m is even and for each k = 1, . . . , d, the number of components of i2 m 
equal to k, denoted by 2mfc, is also even, and a(i n ) = 0, otherwise. 

Lemma 7. We have 

oo 

h = hm{f2m,e)i 
m=0 

where f2m,e is the element of (J-t d )® 2m given by 

f2m,e {hm , r l , ■ • ■ , r 2m) 



(36) 

(2m)! i io ^ 

Proof. Stroock's formula leads to 

f n ,e = ~ t f T f E[Di\';;;Xp £ (Bf - Bf )] ds dt, 
n\ Jo Jo 

where ij 6 {1, . . . ,d}, and r/ 6 [0,T]. 
Let us compute this expectation 

n 

E[D%;;;:£ Pe (B? - Bf)]=E[d^ ■ -eh Pe {Bf - Bf)] ft 1 M (, 

j'=i 

On the other hand, using the Fourier transform 
E[d h ---d in Pe (Bf -Bf)] 

lTT) d jRd 

■J^C 1 --- C n exp{-^(|i - s\ 2H + e)|e| 2 } d£ 



[s,t](r ji 



(27 



i n (2ir)- d / 2 (e + \t- s\ 2H )- d/2 - n/2 E[X h ■ ■ ■ X in ], 
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and the result follows. □ 

Let us compute the L 2 -norm of the nth Wiener chaos of I e . 



E[(/ 2m (/ 2m , £ )) 2 ] = (2m)!||/ 2mi£ ||^ (2m) 

(2m)! (27r) - ( 
(2mi)!---(2m d )! ((2m)! 



(2m). 2^ f . . . Orr, X\ (f0rr,\]\2 "(W 



(37) 



wjiH \-md=iTi 

x / ( e + A)^ /2 - m (e + p)~ d/2 - m /i 2m ^ 



r 

— jf> + X)- d / 2 - m (e + p)- d / 2 -"V 2m dr, 
(2m 1 )!---(2m d )! 



where 



(27r) d 2 2m 

= E 



(mi!) 2 • • • (m^!) 2 



If i7d < 3/2, we can deduce the Wiener chaos expansion of the renormal- 
ized self-intersection local time I: 

oo 

(■ = E ^m(/2m), 
m=l 

where 

m , n , . . . j r 2m J 
(2^/ 2 a(i 2m ) 'T rt 



•/ U J U 7=1 



(2m)! 

5. Proof of Theorem 2. The purpose of this section is to show Theorem 2 . 

Case 3/ (2d) < H < 3/4. We are going to check that the family of ran- 
dom variables e d / 2 ~ 3 /( iH ) (I £ — E(J e )) satisfies the conditions of Theorem 3, 
for any sequence e(k) — > 0. We have 

oo 

£d /2-mm (Je _ m)) = £d /2-m H) ^ hm{hme) . 

m=l 

Thus, from (37) we obtain 

S d ^^M[(hm(f2m,e)) 2 ] 
(2vr) d 2 2m j- 



[(e + A)(e + P )]- d ' 2 ~ m ^ m dT. 
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Making the change of variables yields 

^-3/(2H) E[(/2m(w))2] 



where T £ = e 1 /^ 2H ^T. Using the decomposition T £ n {s < s'} = T{ 

T 3 £ , where % £ = £- l l^% and the sets % are introduced in (46), we obtain 

s d - 3 ^E[(I 2m (f 2m , E )f] 

= , r / (T-(a + b + c)e 1/{2H) ) * m da db dc, 



where 

(38) * m = Et(i + Ai)(i + Pi )r d/2 - m A 2m , 



i=l 



and the functions Aj, p^, /Ltj, have been defined in (47)— (52). 
As a consequence, 

(39) hm^- 3 /^)E[(/ 2m (/ 2m , £ )) 2 ] = (2ff ^Li / R 3 

and this implies condition (ii) of Theorem 3. We also have 

e d - 3/i2H) n(hm(f2 m ,e)) 2 } < (2 ^L-1 Xa *- dQdWC ' 

and, taking into account Remark 1 after Theorem 3, in order to check con- 
ditions (i) and (hi) it suffices to show that 

a 2 := - — r I ^ m dadbdc < oo. 

^ ( 27r )d 2 2m-i y R 3 

From (38) we have 

3 oo „ 

where A j = 1 + Aj , pi = 1 + pi and 
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Hence, 



3 2 f — 

" 2 = g(27F4 (Ai 



Pi 



-d/2 



y (2m)!7f 
^ (m!) 2 2 2 ^ 

. m=0 v ' 



1 



da do dc. 



Using Newton's binomial formula and the definitions (64) and (65) we obtain 

3 



(40) 



Alf (A^)- d/2 [(l -7*r d/2 - l]dad6dc 
(2vr) d ^ -'IRS 

2 3 /" — 
= - — —j / 0j da do dc, 
(27r)^V + 

where the functions 0j are defined in (65). Finally, a 2 < oo by Lemma 13. 

Condition (iv) of Theorem 3 follows from the convergence (43) in Propo- 
sition 8. 



Case H = 3/ (2d). As in the proof of the previous case we need to check 
that the family of random variables (log(l/e)) _1 / 2 (/ e — E(/ e )) satisfies the 
conditions of Theorem 3, for any sequence e(k) — > 0. We have 



(\og(l/e))-^(L - E(J.)) = (log(l/e)r 1/2 £ hM 



2m, e ) 



m=l 



and as before we obtain 
E[(J 2m (/ 2m , £ )) 2 ] 



log(l/e) 



2a r 



(2ir) d 2 2m log(l/e) 



0<a+6+c<Te- 1 /(2ff) 



(T - t(a + b + c)e~ 1/{2H) )^ m dadbdc, 



where ^ m is defined by (38). From the estimates obtained in the proof of 
Lemma 13 it follows that the term (a + b + c)e -1// ( 2 ^ gives no contribution 
to the limit and we can write, using L'Hopital's rule, 

E[(/ 2m (/ 2m>£ )) 2 ] 



lim ■ 

eio log(l/e) 

= lim 



2Ta r 



elO (27T) d 2 2m log(l/e) Jo<a+b+c<Te-VW 



da db dc 



FRACTIONAL SELF-INTERSECTION LOCAL TIME 
2Ta m 



21 



1™ (2vr) a! 2 2m logx J <a+b< v 

v<Tx l/(2H) 



lim 



X l/(2H) T 2 



Too H2 2m {2ir) d Jo <a +b<TxWO 

X 3/(2H) T A 



^ r m (a, b,rj — a — b) dadbd?] 

* m (o, 6, Tx 1/(2H) -a-b)dadb 



lim 



Too H2 2m {2n) d Jo<a+(3< 



* m (Tx 1/(2ff) (a,/3,l - a- (3)) dad/3. 



Set = (a,/3, 1 — a — /?). We have 3/(2H) = d, and 



i=l 

and when x tends to infinity this converges to 

3 



r- 3 x:(Vi)- d/2 7r(^). 



i=l 



Hence, 



(41) 



lim ^[(^ m (/ 2m , £ )) 2 ] 
log(l/e) 



H2 2m (2ir) d ^7o<a+/3<i 

This proves condition (ii) of Theorem 3. 
On the other hand, 



E / (A iPi )- d/2 7r(ea,/3)d«^- 



a" 



3 oo 

-EE 



(\ lPi y d/2 j?(e a ,f3)dadP 



E 



- #(2vr) d 7 <a+/3<l 



E 

. m=0 



(2m)! 7 r 
(m!) 2 2 2 ^ 



(e aj/ g) da d/3. 



Using Newton's binomial formula and the definitions (64) and (65) we obtain 



(J" 



3 

m -uE / (A^)- d / 2 [(l- 7 i)- d / 2 -l](e^)d«d/3 



(42) 



3 

f LE / ((^r d/2 - (X iPi y d / 2 )(e a , p )dadp 
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1 3 f 
= WI^Ad 12 / Qi(e a)/} )dadp. 

Finally, a 2 < oo by Lemma 12. This proves condition (iii). 
Finally, 

-. oo 

OO 3 rri „ 

< E 12 H2 Cr 2ir¥ L ^ amr d/2 T(e a ,p)dad(3, 

m=N+li=l HZ W ^0<a+/3<l 

and this converges to zero as iV tends to infinity. So condition (i) holds. 
Condition (iv) follows from the convergence (44) in Proposition 8. 

Proposition 8. Fix m > 1 and l<p< 2m — 1. 

(i) Ifd>S/(2H), then 

(43) hme 2d ~ 3/H ||/ 2m!e ® p ./WH^^-p) = 0; 

(ii) ifd= 3/(2H) and H < 3/4, tfterc 

(44) lim(log(l/e))~ 2 ||/ 2m , £ ® p /2m, e ||#®2(2 m -p> = 0. 

Proof of (43). The proof will be done in several steps. 

Step 1. Let us first compute f2m,e ®p /2m,e- We have, from (36), 
{hm,e ®p f2m,e)(hm-p,j2m-p,T'l, . . . ,rim~pi w lj • • • ,U2m-p) 

(2tt)^_ 



((2 m )H 2 E a (i2m-p ; kp)a(k p , j 2m - p ) 



jf_(e + \t - s\ 2H )~ d / 2 ~ m {e + \t' - s'\ 2H )~ d / 2 - m ^(s, t, s', t') 



IT 

2m— p 

X II 1 [«,t]( r i) 1 [*' 1 f]( u i) dr - 



As a consequence, 

||/2m,£ ®p /2m,e||jygi2(2m-p) 
4 



(45) =A 2 I f[(e + \t l - Sl \ 2H r d / 2 - m 

x n\ 2 {s, t)nl A (s, t)/J{™~ p (s, t)/i2™~ P (s, *) ds dt, 
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where 



and 



Set 



A2 = ((2m)!) 4 ^ 1 ^ a ( i2m "P' k p)«( k P- J2m- P ) 

s = (s 1 ,s 2 ,s 3 ,S4), t = {ti,t 2 ,t 3 ,U), 
V = {(a, t) G M 8 : < Si < U < T} 

Mi j j ^) — /H^i, ^ij Sjitjj- 



i=i 

x ^ 2 (s, t)/i3 4 (s, t)// 2 ™~ p (s, t)fil™~ p (s, t) ds dt. 
To show (43) it suffices to show that lim £ ^o = 0. 

S'tep 2. We are going to use a representation of <J> £ in terms of the ex- 
pectation of the product of four random variables. Consider a family of 
independent one-dimensional fractional Brownian motions with Hurst pa- 
rameter H , B\ 3 , i = 1, 2, 3, 4, j = 1, . . . , 2m. We can write 



where the product is over the pairs (z,j) = (1,3), (1,4), (2,3) and (3,4), and 

X? = / (e + lt-s^r^U J] (B^-B^)(Bl k -Bl' k )dsdt. 
Jo Jo 1=1 k=i 

Step 3. The variances of the random variables £ d / 2 ~3/( iH ) X l e <i converge 
to a constant times T as e tends to zero. In fact, 



ne d ^ 2H \x^f] = ^^^ 3 /(^)E[(/ 2m (/ 2m , £ )) 2 ] 



and we have already proved in (39) that this converges to 



2T / ^ m dadbdc < oo. 
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Step 4. In order to show that lim £ _»o < ^e = 0, it suffices to show that 
each family of random variables £ d / 2 ~ 3 /( iH ) X % £ converges in law to a normal 
random variable. In fact, taking into account that E(X* , - J X* > J ) = 0, for 
7^ and the results proved in Peccati and Tudor (2003), the four- 

dimensional vector 

e d/2-3/(4fl) (X 1, 3) X 14 X 2 >3) X 3 )4) 

will converge to a Gaussian random vector with independent components 
and, therefore, the expectation of the product of its components will con- 
verge to zero. 

Notice that the families of random variables {e d / 2 ~ 3 /( 4H *> Xpi , e > 0} are 
identically distributed and they have the same distribution as {Z £ ,e > 0}, 
where 

Z£=e d/2-3/(4H) f T [\ £ + \ t - s fHyd/2-rnfr iB H,i_ B H,i )dsd ^ 

Jo Jo i=1 

and {Bf,t > 0} is a (2?n)-dimensional fractional Brownian motion with 
Hurst parameter H. We make the decomposition 

Z e = z\ — z 2 , 



where 

zl = £ mm f (i + \t- s \ 2H )- d / 2 - m z st dsdt, 

zl = £ vm r (i + \t- s \ 2H r d ^ m z st dsdt, 

Ja<tATe-VQH) 

and Z s j has been defined in (16). Let us show that Z 2 tends to zero in L 2 (S7). 
We have 



Z £ 2 || 2 <e 1/(4H) / {l + \t- s\ 2H )- d l 2 - m \\Z st \\ 2 dsdt 

= e V(^) [ (l + \t-s\ 2H y d / 2 - m \t- S \ 2mH dsdt 

Js<Ts-V(. 2H )<t 

< £ v(^) r (i+\t-s\ 2H )- d / 2 d S dt 

Js<Te- 1 /( 2H '><t 

Te -1/(2H) 

= £ V(4H) / / (l + x 2H r d / 2 dxd S 

Jo JTe-VV H )-s 

Te -1/(2H) 

/ / (l + x 2H )- d / 2 dxdy, 

Jo Jv 



which converges to zero because Hd > 3/2. 
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Z. ds ) 



Finally, we have 

,00 / „ N . Te -V(aH3 

Z l = jf (1 + x 2if ) -d/2-m^l/(4H)^ 

which converges in law to a normal distribution from Proposition 5. □ 
Proof of (44). Set 

$ £ = (log(l/e))- 2 f {[(e + lU- Sl \ 2H )- d / 2 ~ m 

As in the proof of (43), it suffices to show that <3> £ tends to zero as e tends 
to zero. We can also write 



$ £ =(iog(i/ e ))- 2 E^n^'), 



where the Xp J have been defined in step 2 of the proof of (43). 

From (41) it follows that the variances of the random variables (log(l/e)) -1 / 2 
X^' 3 converge to a constant. In fact, 



E[(log(l/ e ))-W) 2 ] = LJ (log(l/e))- 1 E[(/ 2m (/ 2m , £ )) 2 ]. 



Then, as in the proof of part (i), it suffices to show that the family of 
random variables (log(l/e)) _1 / 2 X*'- ? converges in law to a normal random 
variable as e tends to zero. This family has the same distribution as {Y e ,e > 
0}, where 



pT pt 2m 

y £ = (log(l/e))- 1 /2 / / ( e + \t-s\™)-^"I[{B*> i -Bf< i )d8dt, 
Jo Jo f- = \ 

and {Bf*,t > 0} is a (2m)-dimensional fractional Brownian motion with 
Hurst parameter H. As in the proof of part (i) we make the decomposition 



Y —Y 1 —Y 2 



where 



y £ 2 = (loga/e))- 1 / 2 ^ 4 ") / (i + \t- s \ 2H r d / 2 - m z st dsdt, 

Y} = (log(l/e))- 1 / 2 e 1 /( 4 ^ / (1 + \t-s\ 2H )- d ' 2 - m Z st dsdt, 
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and Z s j has been defined in (16). As in the proof of part (i) we obtain 

Te -1/(2H) 

||y £ 2 || 2 <(log(l/ £ ))- 1 /2 £ V(4H) / / (l + x 2 ^)-^ 2 ^, 

JO Jy 



which clearly converges to zero. Finally, the family 
= (logU/e))- 1 / 2 J^il+x™)-^-™ ^ 1 ^ f 



Te -1/(2H) 



Z 



s,s+x ' 



ds dx 



converges to a normal distribution from Proposition 6. □ 



Remark 3. The limit variance appearing in (40) can be also expressed 



as 



a' 



(27T) 



l + x 2H )(l+y 



2H\ 



{\z + y\ 2H + \z-x 



2H 



-d/2 



{{1 + x 2H ){1+y 2 H)) - d /2 



dx dy dz. 



We do not know a simple expression for this variance, even in the case 
H = 1/2. 

Remark 4. In the particular case H = 1/2, and d = 3, the limit variance 
appearing in (42) is a 2 = (27T 2 )" 1 . In fact, in this case we obtain 



a 



1 

47T 3 



0<a+6<l 



[(ab+(a + b)(l-a-b)y 3/2 

((a + 6)(1 - a))" 3/2 + (6(1 - 6))" 3/2 - b^' 2 } dadb. 



Making the change of variables x = a/(a+b), y = l/{a + b), it is not difficult 
to check that the first integral equals 2ir. The remaining terms are equal to 
Jo 6 -3 / 2 (l + b — (1 — b)' 1 / 2 ) db, and with the change of variable \J\ — b = 
cos#, it is clear that this integral vanishes. 



6. Technical lemmas. We will denote by k a generic constant which may 
be different from one formula to another one. We will decompose the region 
T defined in (9) as follows: 



(46) 



T n{s<s'} = T 1 UT 2 UT 3 , 
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where 



T 1 = {(t,s,t', S ') 
T 2 = {{t,s,t',s') 
T 3 = {(t,s,t',s') 



0<s<s' <t<t' <T}, 
0<s<s' <t' <t<T}, 
0<s<t<s' <t' <T}. 



We will make use of the following notation: 

(i) If (t, s, t' , s') 6 T±, we put a = s' — s, b = t — s' and c = t' — t. On this 
region, the functions A, p and p defined in (10) and (11) take the following 
values: 



\2H 



(b + c) 
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u2H C 2H 



-.2Hi 



(47) \ = \i:=\ 1 (a,b,c):=(a + by~, p = p x 

(48) fi = p! :=n 1 (a,b,c) := \[{a + b + c) 2H + b 2 

(ii) If (t, s, t', s') £T 2 , we put a = s' — s, b = t' — s' and c = t — t'. On this 
region we will have 

p = p 2 :=(a + b + cf H , 



(49) 
(50) 



A = A 2 := b z 



p = p 2 :=\[{b + c) 2H + {a + b) 



\2H 



„2H 



2H] 



(iii) If (t, s, t' , s') £ 73, we put a = t — s, b = s' — t and c = t' — s' . On this 
region we will have 



2H 



„2H 



(51) A = A 3 := a"' , p = p 3 := c ~ 

(52) ti = fX3 :=^ a + b + c fH + b 2H_ {b + c) 2H , ^ 



For 1 = 1,2,3 we set 5, = Aj/?j — p? Oj = 5. 



-d/2 



(a + 6)^ 

(AiPi)"^ 2 . Note that Aj, 
Pi, /ij and so on, i = 1,2,3, are functions of a, b and c. 

The following lower bounds for the determinant of the covariance matrix 
of B?' 1 - Bf' 1 and B*?' 1 - Bf l , were obtained by Hu (2001) using the local 
nondeterminism property of the fractional Brownian motion [see Berman 
(1973)]. 

Lemma 9. There exists a constant k such that, for all a,b,c> 0: 
(i) 

(53) 6x > k[(a + b) 2H c 2H + {b + c) 2H a 2H ]. 
(ii) Fori = 2,3 

(54) 5i > kX iPi . 



The following lemma provides some useful inequalities. 
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Lemma 10. We have 

(55) k(a + b + c) 2H ~ 2 ac< P3 <kb 2H ~ 2 ac. 
For i = 2, 3 we have 

(56) e, < kn 2 {\ iPi )~ d i 2 - 1 

and 

(57) e, < k(\ lPt y d / 2 . 

Proof. The inequalities in (55) follow from 

^ = \ [(a + b + c) 2H + b 2H - (a + b) 2H - (b + c) 2H ] 

= H{2H-l)ac f 1 f\b + vc + ua) 2H ~ 2 dudv. 
Jo Jo 

We have, for i = 2, 3, 

r / ..2 \ -d/2 n 

(A^r d/2 . 



2 x -d/2 
1-^-1 "I 



The estimate (54), assuming k < 1, implies j^j- < 1 — A; and (56) holds. 
Moreover, (54) also implies (57). □ 

The following lemma is the basic ingredient in the proof of Theorem 1. 

Lemma 11. Let be defined by (15). Then Et < oo if and only if 
dH < 3/2. 

Proof. 

Step 1. Suppose dH < 3/2. We claim that 



(58) / @i da db dc < oo 

J[0,T]3 

for i = 1,2, 3. From (53) we deduce 

8 1 >k(a + b) H (b + c) H a H c H 

(59) 

>k(abc) iH '\ 
Then, (59) together with the estimate 

\ lPl = ( a + b) 2H {b + c) 2H > (abc) 4H / 3 
implies (58) for i = 1. 
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To handle the case % = 2, we decompose the integral over the regions 
{b > ija}, {b > rjc} and {b < rja, b < r/c}, for some fixed but arbitrary rj > 0. 
We have, using (57), 



O2 dadbdc < k 



da db dc 



lb> V a - Jb> V a(a + b + c) dH b dH ' 

If dH < 1, then this integral is finite. If 1 < dH, then 

- Jo Jo (a + c) dH Jr,a 

<k( I a- 4dH / 3+1 c- 2dH ^dadc < 00. 



JO 

It is also easy to show that S < 00 in the case 1 = (iff. The case b>r/c can 
be treated in a similar way. 

To deal with the case both b < rja and b < rjc, we make use of the esti- 
mate (56) and the following upper bound for (j^: 

A*2 = o K a + 0) ~ a + (0 + c) - c J 

(60) 

for 77 small enough. In this way we obtain 

6 2 < k(a AH ~ 2 + c AH - 2 )(a + b + c) -^-dH b 2-2H-dH 

< k[a^ d ^ H b dH / 3 + c (2-d/3)^ 6 «w/3] ( a + 6 + c) -2H-ciff 6 -dff_ 

Hence, 



Jb<ria,t 



©2 <ia c?6 dc 

' b<rja,b<Tjc 



-dH/ , 1 , „\-2H-dH 



< k / 6^" (a + 6 + c) 

J b<r]a,b<rjc 

<jfc / 6" dH (a + 6 + c)^ 2H " dH a( 2 - d / 3 ) H ^/ 3 (ia(i6d C 

J[0,T]3 

< it / b -2dH/3 a -2dH/3 c -2dH/3 da db ^ 

J[0,T]3 

which is finite if <ii7 < 3/2. 

To handle the case i = 3, we decompose the integral over the regions 
{a>rjib,c>f]2b}, {a < 771 6, c < {a > 7716, c < 772^} and {a < 7716, c> 
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772 fr}. By symmetry it suffices to consider the first three regions. We have, 
using (57), 



G3 dadbdc< k db 

> JO Jn 



T da f T dc 



a>r]\b ,c>r)2b JO Jr)\b Jr)2b C^ 1 ^ 

Let us now suppose both a < rub and c < r/2&. Using (55) and (56) and that 
H < I yields 

G 3 < kb m - i a 2 - 2H - dH C 2 - 2H ~ dH < ka -2dH/3 c -2dH/3 b -2dH/3^ 

which implies that the integral over this region is finite. Finally, let us con- 
sider the case c<rjib and a > r\ib. If Hd > 1, then (57) yields 

-Hd 



63 < k{acY 

<r _ 

d — 2' 

^ = 1 [ (a + b + c) 2H - (a + b) 2H - (c + b) 2H + b 2H ) 



which is integrable. So, we can assume H < 4 < 4. Then 



(61) 

<kb 2H ~\ 

if 772 is small enough. Hence, using (56) and (61) we get 

e 3 < kb 4H ~ 2 a~ 2H - dH c 2 - 2H - dH . 

Consequently, if -dH + 2H + 1 < 0, 



/ 



O3 da db dc 

I c<r]ib,a>T]2b 

b AH-2 a -2H-dH c 2-2H-dH dcdbda 



c<riib,a>r]2b 

<k [ a- dH ~ 2H b- dH+2H+1 dbda 

<k [ T a~ 2dH+2 da, 
Jo 



which is finite if dH < 3/2. The case —dH + 2H + 1 > is easier. 

Step 2. Assume that H = 3/ (2d), and let us show that = 00. It suf- 
fices to prove that 

(62) ~ — ~ —-7 / ii 2 (X P )- d l 2 - l dsdtds'dt' = 00, 



2(2vr) d Jr 
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because from the identity (37) this is the second moment of the second chaos 
of the renormalized self-intersection local time. In order to check (62) we will 
show that 

A:= f {T-a-b-c)nl(X 3 p 3 y d/2 ~ 1 dadbdc = oo. 

Jo<a+b+c<T 

With the above notation, we have, using (55), for e > small enough, 
A>k ( {b + c + a) iH - 4 {ac) 2 ~ Hd - 2H dadbdc:=B. 



If d = 2, we get 



B = k - — dadbdc = oo. 

'[o,e] 3 (a + 6 + c)ac 



For d > 2, we have 2 - Hd - 2H = 1/2 - 3/d > -1. Hence, 

B = ^JT7 I l(c + a) 4H - 3 -{e + c + af H -\acf- Hd ~ 2H dadc, 
3 - AH J[o, £ ]2 



and 



f {c + af H -\ac) 2 - Hd - 2H dadc 

J0<a<c<e 

> 2^H-3 f „2-Hd-2H JZH-l-H. 



0<a<c<e 

> k I" a 2 ~ 2Hd da = oo 



because 2H — 1 - Hd < - 1 . □ 

Lemma 12. If H = 3/ (2d) < 3/4, then for i = 1,2, 3, 

/ / &i(a, b — a, 1 — b) dadb < oo. 
Jo Jo 



Proof. Suppose first i = 1. Then, using (53), we obtain 

/ / 5 1 d ^ 2 (a,b — a,l — b) dadb 
Jo Jo 

<k f 1 t '[6(1 - b) + (1 -a)a]" 3/2 dadb < oo 
Jo Jo 



and 



1 r b 

(X lPl )- d/2 {a,b - o,l - b) dadb 



o Jo 



fl 

< k / b~ 3/2 [1 - (1 - b)~ 1/2 } db < oo. 
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For i = 2 we can write, using (56), 



n©2(o> b — a, 1 — b) dadb 



<k f 1 (\b-a)-*/ 2 ~ 2H [{l-a) 2H + b 2H -{l-b) 2H -a 2H ] 2 dadb 
Jo Jo 

<k f 1 [ b (b-a)-V 2 - 2H (b 2H -a 2H ) 2 dadb. 



10 Jo 

Then, if rj is small enough, 



f ( p_ a) -3/2-2H( b 2H_ a 2H ) 2 dadb 
J a<b,b—a<rja 

<k f (b - a) 1 ' 2 -™ 'a m ~ 2 dadb < oo 

J a<b,b—a<r]a 

because H < 3/4, and 

f {b _ a) -3/2-2H {b 2H _ a 2 H) 2 dadb 

J a<b,b—a>na 



' a<b,b—a>r]a 

<k ^ x-V 2 - 2H+2MH dadx 



JO 

k C x- l ' 2 - 2H+2MH dx<oo. 



Finally, for i = 3, we obtain, using (56), 



n03(a, b — a, 1 — b) dadb 



<k f 1 [ b a- 3 / 2 - 2H (l-b)- 3 / 2 - 2H 
Jo Jo 

x [1 + (b - a) 2H - (1 - a) 2H - b 2H } 2 da db < oo. 

In fact, if we fix 771 and 772 small enough, it suffices to check that the above 
integral is finite over the regions 

G\ = {a < b, a > f]i(b — a) , I — b > r]2(b — a)}, 

G2 = {a < b, a < 771(6 - a), 1 - b < 772(6 - a)}, 

G 3 = {a <b,a> 771(6 - a), 1 - b < 7/2(6 - a)}, 

G4 = {a < b, a < 771(6 — a), 1 — 6 > 772(6 — a)}. 

On C?i we have, using (57), 

/ @ 3 (a,b-a,l-b)dadb<k f a~ l / 2 ~ 2H [l -b)- 1,2 ~ 2H dadb. 
Jd 7gi 
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In order to estimate this integral we make the decomposition 

Gi C (G?i n {a > e}) U (Gi n {1 - b > e}) U (G x n {a < e, 1 - 6 < e}). 

On the set Gi H {a > e} the integrand is bounded by fc(l — b)~ l l 2 ~ 2H , and 
making the change of variable y = l — b, x = b — a, we obtain a finite integral 
because H<\. The same approach is used to handle the integral over G± n 
{l-b>e}. Finally, 

Gi n {a < e, 1 - 6 < e} C Gi n {fe - a > 1 - 2e}, 

and, on this set the integrand is bounded. 
From (55), we obtain on G2 

1 + (ft _ a fH _ (1 _ fl) 2H _ b 2H < fc(fe _ a) 2ff-2 a(1 _ 6) 

Hence, 

/ e 3 (a,b-a,l-b)dadb<k f a 1 ' 2 -™ {I - b) 1/2 ~ 2H {b - a^^dadb. 
J G2 JG2 

On the set G2 H {6 — a > e} this integral is clearly finite, and on the set 
G2 H {b — a < e} we have b — a > 1 — 2s, and again the integral is finite 
provided e < 1/2. 

Fix < £ < 1 — 7/2. Using (56), on the set G3H {b — a > e} we have 

e 3 ( a , 6 _ a , 1 - 6) < ka x l 2 - 2H (1 - 6)V2-2H £ 4H-4 

which is integrable. On the set G3 n {6 — a < e} we have 

a = 1 - (6 - a) - (1 - b) > 1 - e - r/ 2 > 0. 

Then, using the estimate 

1 + ( 6 _ a )2" _ (1 _ a fH _ h 2H<\ 1 _ 6 2H| + | (6 _ a) 2H _ (1 _ a) 2H| 

we obtain, on the set G3 n {b — a < e}, 

e 3 (a, 6 - a, 1 - b) < ka~ 3 / 2 - 2H (l - b)^ 2 - 2H [(b - af H ~ 2 + a m ~ 2 } 
<k(l-b) 1 / 2 - 2H [(b-a) 4H - 2 + l), 
which is again integrable. □ 

We will make use of the following notation, for i = 1,2, 3: 

(63) Xi = Xi + 1, p i= Pi + 1, 

(64) 5i = XiPi - fj, 2 , 

(65) e i = (ff i )-* /2 -(Ai&r ,i/2 - 
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Lemma 13. Letd>2 and 3/ (2d) < H < 3/4. Then for i = 1,2, 3 
(66) / Qidadbdc<oo. 



Proof. The proof will be done in several steps. 

Step 1. Suppose first i = l. Then, 

0! < {hy d/2 = (<Si + (a + 6) 2 " + (6 + c) 2 " + l)~ d/2 . 

First, Gi is integrable in the region {a, b, c > 1} because the estimate (59) 
yields 

G\ < 5 1 " d/2 < fc(a&c)- 2<Hr / 3 
and dH > 3/2. On the other hand, using (53) we obtain 

@i < [k(a + b) 2H c 2H + k(b + c) 2H a 2H + (a + 6) 2H + (6 + c) 2H + l]~ d / 2 
< fc[max(a 2H + b 2H , a 2H + c 2 ", b 2H + c 2H ) + l]" d / 2 , 
which implies the integrability of Qi on the complement of the region {a, b, c > 

S'tep 2. Suppose that i = 2. Then, using (54) we obtain 

Q 2 <{5 2 + {a + b + c) 2H + b 2H + l)- d ' 2 

<[k( a + b + c) 2H b 2H + (a + b + c) 2H + b 2H + l]~ d / 2 . 

On {a, b, c < 1} the function O2 is clearly integrable. If one of the coordinates 
is bigger than 1 and the other two are smaller we use the estimate 

Q 2 < A;[max(a, b,c) 2H + l]" d/2 . 

If a < 1, b > 1 and c > 1, or if a > 1, b > 1 and c < 1 we use the estimate 

Q 2 < A;[max(a6,6c) 2// + l]- d / 2 . 

If a > 1 , b > 1 and c>l we use the estimate 

§ 2 < k[{a + b + c) 2H b 2 H] -d/2 < fc ( a6c )-W/3_ 

So it only remains to check the integrability of G2 in the region 

G = {(a,b,c) :b < l,a,c> 1}. 
Using inequality (56), on the region G we have 
G 2 < kn 2 (\ 2 p 2 )- d ' 2 

< fc [(l + c )2^ + (o + !)2H _ C 2H _ a 2 H] 2 {a + c yH d -2H 
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Consequently, 

O2 da db dc 



G 



< kj™ y°°[(a + l) 2H - a 2H ] 2 (a + c y™~™ dadc 

< k J °° J °° a 4H ~ 2 (a + c) ~ Hd ~ 2H da dc 
<kj™ a 2H - Hd ~ 1 da<oo 



because d > 2. 



Step 3. Suppose that i = 3. We claim that (66) holds for i = 3. On 
{a > 7/1 6, c > r?2&} we have using (54) 

e 3 <(52 + a 2H + c 2H + l)- d / 2 

<k(a 2H + l)- d ' 2 {c 2H + l)- d ' 2 1 

and as a consequence, we obtain 

03 da db dc 

a>r)ib,c>r)2b 

db Ll lb a ~ Hdc ~ Hddadc - k k 6^ <CX) ' 

if dH>3/2. 

On {a < r/16, c < 7/26}, using the estimates (56) and (55), we obtain 
©3 < kb iH -\ac) 2 {a 2H + iy d l 2 -\c 2H + l)-^ 2 - 1 . 
If dH > 3/2 and H < 3/4, then 

©3 da db dc 

a<r]\b,c<r]2b 

<k f (a 2H + iy d l 2 - l (c 2H + l)- d l 2 - l (ac) 2 dadc [°° b m ~ A db 
J[0,oo) 2 Ja/r)iVc/iK 

<k[ (a 2H + iy d /^\c 2H + ir d / 2 ' 1 (ac) 2 c 4H - 3 dadc<oc. 

Ja<c 

Finally, let us consider the case c < 7/26 and a > r/16. We have, using the 
estimates (61) and (56), 

©3 < Kb 4H ~ 2 + (a + b) AH - 2 )c 2 (a 2H + l)- d ' 2 - l (c 2H + 



3G 



Y. HU AND D. NUALART 



if rji is small enough. If 2H < 1 we use (a + b) iH 2 < b m 2 , and if 2H > 1 
we use 

(a + 6) 4ir-2< (1 + ??2) 4if-2 6 4H-2 

In this way we obtain 

6 3 < k{b AH ~ 2 + a 4H ~ 2 )c 2 (a 2H + ir d / 2 -\c 2H + l)-^ 2 - 1 . 

By decomposing the integral over the regions {1 < c < r]2b < ^a} , {c<l< 

V2b < ^" a } 5 { c < V2b < 1 < ^-a} and {c < 7726 < < 1} we easily show that 
the integral is finite. □ 
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